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ABSTRACT 


Rate  coefficients,.,**1*  calculated  for  the  reactions  of 
hydrogen,  deuterium,  and  deuterium  hydride  with  chlorine. 
The  potential  energy  surfaces  of  the  isotopic  transition 

'  i  jS  *  ✓  ' 

complexes  of  HCl^  and  H^Cl  have  Won  constructed  bythe 
Sato  method  to  match  experimental  activation  energies  in 
accordance  with  absolute  rate  theory.  The  rate  coeffi¬ 
cients  calculated  using  absolute  rate  theory  are  then  used 
to  evaluate  the  relative  rates  of  reaction  of  H^/HD,  H^/D^ 
with  chlorine  atoms,  and  the  relative  rates  of  the  fast  re¬ 
action  step  H  +  C  1  ^  compared  to  the  back  reaction  H  +  HC  1 . 
The  results  agree  very  well  with  the  few  experimental 
values  reported  in  the  literature.  Otrr  theoretical  calcu¬ 
lations  show  that  prior  experimental  values  of  the  pre¬ 
exponential  factors  in  the  Arrhenius  -  type  expression  for 
rate  constants  are  not  required  in  order  to  obtain  better 
information  on  the  structure  and  vibrational  frequencies  of 
the  transition  states.  The  agreement  of  the  activated  com¬ 
plex  theory  with  experiment  was  Tound  be  satisfactory. 
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INTRODUCTION 

Jones  *  measured  the  relative  rates  of  reaction  of  H  and  HT  with 

2  ^ 

chlorine  atoms,  and  Ashmore  and  Chanmugan  measured  the  slow  reaction 
between  hydrogen  molecules  and  chlorine  atoms.  Measurements  have  been 

made  of  the  relative  rates  of  reaction  of  chlorine  atoms  with  H-,  and  D,  by 

3  C  4  L 

Rollefson  ,  with  mixtures  of  H7,  HD,  and  D?  by  Bigeleisen,  et  al,  and 

5  ^  h 

Farkas  and  of  Cl^  with  H,  and  H  with  HC1.  In  Fig.  1  some  of  these  experi¬ 
mental  data  are  compared  with  values  calculated  by  Wheeler,  Topley,  and 
Eyring"^  using  the  London  potential.  The  value  of  the  largest  discrepancy  in 
the  ratio  of  the  experimental  relative  rates  to  these  theoretical  relative  rates 
is  about  five.  In  addition,  the  difference  in  theoretical  activation  energies 
for  the  reactions  of  H^  with  Cl,  and  HT  with  Cl,  has  the  opposite  sign  from 
that  of  the  experimental  data. 

g 

Bigeleisen  and  Wolfsberg  obtained  satisfactory  agreement  between 

theory  and  experiment  by  using  the  experimental  value  of  the  pre- exponential 

factor  for  the  relative  rates  of  reaction  of  mixtures  of  H^,  HD,  and  D^  with 

Cl  to  estimate  the  symmetrical  stretching  vibration  in  the  H_C1  transition 

6  ^ 

state.  More  recently,  Klein  and  Wolfsberg  have  used  a  similar  procedure 
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to  estimate  the  doubly  degenerate  bending  vibration  in  the  HCl^  transition 
state.  These  estimated  frequencies  are  based  on  low  temperature 
measurements . 

In  this  study,  the  fundamental  vibrations  of  each  activated  complex 

are  calculated  without  prior  information  on  the  experimental  values  of  the 

pre-exponential  factors.  The  molecular  data  for  the  activated  complex 

9 

obtained  from  the  potential  energy  surfaces  using  the  Sato  method  ,  along 
with  molecular  structure  data  for  the  reactants,  are  used  to  calculate  the 
absolute  reaction  rates  of  hydrogen,  deuterium,  and  deuterium  hydride  with 
chlorine.  The  theoretical  results  are  compared  with  experimental  values 
for  the  relative  rates  of  reaction  of  H^/HD,  H^/D^  with  chlorine  atoms,  and 
for  the  rate  of  the  fast  reaction  step  H  +  Cl^  relative  to  that  of  the  back 
reaction  H  +  HC1. 

COMPUTATIONAL  PROCEDURE 
Consider  the  elementary  exchange  reaction 

kf 

A  +  BC  J  AB  +  C  (1) 

The  theory  of  absolute  rates  ^  gives  the  following  formula  for  the  rate 
coefficient  of  this  elementary  exchange  reaction: 

kf  =  *T(?)o^exp(-Eo/kT)  (2) 

In  Eq.  (2)  the  subscripts  c,  a,  and  m  refer  to  the  activated  complex,  the 
reactant  atom,  and  the  reactant  molecule,  respectively;  k  is  the  transmission 
coefficient,  Q  is  the  total  partition  function,  r  is  the  tunnel  effect  correction 
factor,  Eq  is  the  classical  activation  energy  with  zero-level  correction,  and 
the  other  symbols  have  common  meanings.  Substitution  of  the  proper 
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partition  functions**  for  both  complex  and  reactants  allows  one  to  write 
Eq.  (2)  as 


kf  =  A(T)exp[-E0/(kT)] 


(3) 


where  the  pre-exponential  factor 


Am  =  "(s^)  ('  -  «xp[-(hv)/(kT)])  ,n  (l  -  .xpHhVjt/kT) 

and  m  is  the  mass,  I  is  the  moment  of  inertia,  <r  is  the  symmetry  number, 
v  is  the  vibrational  frequency  of  the  diatomic  molecule,  and  is  the  fre¬ 
quency  of  the  complex.  The  moment  of  inertia  of  the  complex  and  the 
vibrational  frequencies  of  the  complex  are  calculated  from  the  potential 
energy  surface  constructed  to  interpret  the  chemical  reaction.  One  must 
assume  a  priori  some  structure  for  the  transition  complex  before  the  potential 
energy  surface  can  be  constructed.  For  the  computations  made  in  this  study, 
a  linear  structure  is  assumed  for  the  activated  complex  since  the  smallest 
increase  in  the  potential  energy  of  the  system  results  if  all  three  atoms 
remain  on  a  straight  line.  The  vibrational  frequencies  of  the  transition 
complex  are  obtained  from  the  force  constants  of  the  energy  surface  by  the 
method  of  small  vibration  analysis.  For  a  linear  transition  state  complex 
consisting  of  three  atoms,  there  are  one  symmetric  stretching  frequency 
vg,  one  bending  frequency  which  is  doubly  degenerate,  and  an  imaginary 
assymetric  stretch  frequency  v^.  The  imaginary  vibrational  frequency 
affects  the  reaction  rate  only  through  the  tunnelling  correction  factor  t.  The 
tunnelling  correction  is  quite  large  at  room  temperature  but  is  less  signifi¬ 
cant  at  temperatures  in  excess  of  500°K. 
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Recently,  Sato  proposed  a  new  method  for  drawing  the  potential 
energy  surface  required  to  interpret  chemical  reactions.  This  method  is 
not  an  absolute  one,  since  one  parameter  K  (the  square  of  the  overlap 
integral)  is  an  adjustable  parameter.  The  Sato  potential  gives  results 
slightly  different  from  the  London  potential*2  used  by  Eyring  and  Polanyi. 

The  essential  differences  between  the  Sato  method  and  the  Eyring  method  of 
constructing  potential  energy  surfaces  have  been  discussed  adequately  in  a 
recent  paper  by  Weston*  . 

An  electronic  computer  was  used  to  calculate  the  potential  energy 
surface  of  a  system  of  three  atoms  using  the  Sato  method  of  constructing 
potential  energy  surfaces.  In  a  system  consisting  of  three  atoms,  if  one  atom 
is  kept  at  infinite  separation  from  the  other  two,  then  two  equations  for  the 
energy  of  the  bonding  and  antibonding  states  derived  by  the  Heitle r- London 
method  are: 

(for  bonding  states) 

Eb  =  (Q..  +  c..)/(l  +  S2.)  =  (Q..  +  c..)/(l  +  K)  (5) 

lj  i)  ij  '  ij'  ij  i) 

(for  antibonding  states) 

E*j  =  (Qij  -  -  sfj)  =  <°y  -  ■„>'<>  •  K>  (6) 


In  these  equations  Q^.  is  the  Coulomb  energy  for  atoms  i  and  j,  a.j  and  S.j  are 
the  corresponding  exchange  integral  and  overlap  integral,  respectively,  and  K 
is  the  square  of  the  overlap  integral.  The  energy  lor  the  bonding  state  was 
calculated  by  using  the  More  potential  function 


Ey 1  ■  ry)]  - 2  HrM'u  ■  ru>] 


(7) 


where  D. .  is  the  dissociation  energy,  r..  is  the  distance  between  atoms  i  and 
0  13 

j,  r„  is  this  distance  in  the  equilibrium  configuration,  and  p„  is  a  constant 
obtained  from  spectroscopic  data.  For  the  energy  of  the  antibonding  states, 
Sato  proposed: 


lexp[-2Pij(rij  -  r?jfl +  2  exp[-pij(rij  -  rSj)]| 


This  relationship  for  the  energy  of  an  antibonding  state  agrees  well  with  the 

calculated  data  for  the  energy  of  the  repulsive  state  of  hydrogen  obtained  by 
14  15 

more  exact  methods  ’  .  If  one  assumes  K  is  a  constant  for  a  specified 

system,  then  Eqs.  (5)-(8)  can  be  solved  for  Q..  and  a.,  as  a  function  of  r... 

XJ  XJ  12  XJ 

The  total  energy  for  each  configuration  is  obtained  from  the  London  equa¬ 
tion  for  a  system  of  three  atoms 


E  =  ( I  +  K)  Qi2  +  Q13  +  Q23  *  2  (°12  “  a13>  +2(a12'a23) 


,  1.  >2  2 
+  2  °13  "  a23 


The  saddle  point  locus  and  three  force  constants  fj,  f^,  and  f^  corre¬ 
sponding  to  the  activated  complex  were  evaluated  by  the  method  of  least 
squares.  The  bending  force  constant  f  by  definition  is 


■(4) 

\90  / r,,,r 


12  23 
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Since  and  0  are  related  by  the  cosine  law,  it 
bending  force  constant  for  a  linear  complex  is: 


is  easy  to  show  that  the 


f»  =  _  _1_  /9E\ 
r12r23  r13  \a02/r12, r23 


(11) 


and  by  use  of  Eq.  (9) 


.12) 


By  the  method  of  small  vibrations^,  the  frequencies  for  a  transition 
state  composed  of  three  atoms  become: 


_  — 
s  2ttc 


_ 

^b  2ttc 


_ 

Vu  2irc 


q  +  ja2  -  4MjM2M3M(f1f2  -  f12)  * 


2MjM2M3 


-.1 

2 


/ — 2 — \ 

Vl2r23/  r12r23  '  1 


a  -  |a2  -  4M1M2M3M(fJf2  -  f^jj2 


2M  j  M2M3 


(13) 


(14) 


(15) 


i 
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t 

where 


a  =  Mj(M2  +  M3)£j  +  M3(M2  +  Mj)f2  -  2MjM3f12 


The  moments  of  inertia  of  the  complex  and  for  the  molecule  are  given  by 


I 

c 


Mj(M2  + 


M3)ri2  + 


2MlM2r12r23 


m3(mt 


M2)r23. 


(16) 


T  M2M3  /  0  \2 
m  "  M2  +  M3  \  2  3/ 


(17) 


In  one  dimension,  the  energy  surface  of  the  reaction  path  has  a  negative 
curvature  corresponding  to  a  maximum  potential  energy.  For  reactions 
involving  light  particles,  penetration  through  the  energy  barrier  is  significant 
below  500°K.  The  quantum  correction  for  penetration  of  such  particles 
through  the  energy  barrier  is  called  the  tunnel  effect.  Using  the  WBK 
approximation,  Bell  * 7  has  recently  obtained  the  permeability  G(W)  of  a 
parabolic  potential  barrier  as  a  function  of  the  energy  W  of  the  incident 
particle  of  mass  m: 


G(W)  =  [1  +  exp(py)] 


(18) 


where 


y  =  i 


w 

E 


(3  =  2ir^a(2mE)  *  ^/h 
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and  E  and  2a  are  the  height  and  width  of  the  barrier.  The  quantum  correction 
for  the  reaction  velocity  is  obtained  by  averaging  the  permeability  G(W)  over 
Maxwell -Boltzmann  energy  distribution 


=  exp(E/RT)  f  exp(-  W/RT)G(W)dW 


_n_o/P 

sin  if  a/p 


a  exp(c  -  p) 


CD 

♦£ 

n=  1 


(-1) 


(n  + 


P  -  a 


exp  (-nP) 


(19) 


where  a  =  E/kT. 

The  permeability  G(W)  is  a  poor  approximation  when  W  ~  E  because  the 
WBK  approximation  is  valid  only  when  W  «  E.  This  approximation  is 
therefore  valid  only  when  the  degree  of  tunnelling  is  large.  The  assumption 
that  G(W)  =  1  for  W  >  E  neglects  the  reflection  of  particles  which  takes  place 
when  the  value  of  W  ^  E. 

If  one  writes 


u  =  2WP  =  -  ihvu/kT  (20) 

then  the  Bell  solution  for  the  tunnelling  correction  effect  reduces  to  the 
Wigner  u  / 24  correction  for  small  tunnelling.  The  tunnelling  correction 
factor  is  evaluated  by  calculating  the  value  of  E  corresponding  to  the  best 
parabolic  fit  to  the  Sato  energy  barrier  over  which  there  is  appreciable 
tunnelling.  This  procedure  is  a  better  approximation  to  the  parabolic  barrier 
than  is  the  parabolic  fit  obtained  from  the  harmonic  frequency  corresponding 
to  the  curvature  at  the  top  of  the  potential  energy  barrier.  At  temperatures 
below  500°K,  this  tunnelling  correction  is  usually  very  large  and  most  calcu¬ 
lations  indicate  that  the  potential  barrier  calculated  by  the  Sato  method  is  too 
narrow.  With  such  large  values  of  the  tunnel  correction,  one  finds  a 
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non-Arrhenius  dependence  in  both  the  absolute  and  relative  rates.  This  does 
not  agree  with  the  experimental  facts.  Subsequently,  we  calculated  rate 
constants  which  included  the  tunnel  effect  at  temperatures  above  500°K  and 
extrapolated  for  the  values  of  the  rate  constants  below  500°K. 

RESULTS  AND  DISCUSSION 

The  following  elementary  reactions  can  occur  between  chlorine  and  a 
mixture  of  H^,  HD,  and  D^: 


2 1 

H  +  Cl  HC1  +  H 

L  k 

K2 1 


22 

H  +  Cl,  Tt  HCl  +  Cl 
2 

k22 


23 

D,  +  Cl  ~ *  DC1  +  D 
k23 


24 

D  +  Cl,  72  DC1  +  C1 
k24 


25 

HD  +  Cl  72  H  +  DC  1 
*25 


*26 

DH  +  C 1  2^  D  +  HC 1 
*26 


(21) 


(22) 


(23) 


(24) 


(25) 


(26) 
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1 9  .  « 

The  isotopic  competitive  method  is  usually  used  to  measure  the 

relative  rates  of  isotopic  reactions.  In  the  case  of  reactions  (21)  and  (22)  it 

is  very  difficult  to  reproduce  the  chain  length  experimentally.  Unless 

reaction  (22)  can  be  quenched,  an  experiment  designed  to  measure  the 

relative  rates  of  reaction  of  H^/HD  with  Cl  atoms  will  in  reality  measure 

R  =  ^21^^25  +  ^26^'  Although  measurements  of  the  individual  ratios 

k^,/koc  and  k,  ,/k-,  would  be  preferable,  this  is  difficult.  However, 

C  1  CD  C  1  4  co 

Bigeleisen,  et  al  ,  evaluated  the  relative  rate  R  by  intercomparing  the 
isotopic  composition  of  the  unreacted  hydrogen  at  measured  amounts  of 
reaction.  A  correction  was  made  by  Bigeleisen,  et  al,  for  the  presence  of 
which  could  be  formed  by  the  backward  reactions  (23)  and  (24). 


Portions  of  the  Sato  potential  energy  surfaces  for  reactions  (21)-(26) 
are  shown  in  Figs.  2-7.  These  potential  surfaces  are  for  linear  configu¬ 
rations  A-B-C  with  r^,  r^^  as  the  coordinates  between  atoms  A  and  B,  and 
B  and  C,  respectively.  In  these  plots  the  contour  lines  are  expressed  in 
units  of  kcal/mole.  The  experimental  activation  energy  of  reactions  (21)  and 
(22)  were  used  to  determine  the  parameter  K.  The  same  value  of  K  was  used 
for  isotopic  modifications  of  reactions  (21)  and  (22). 


The  pre- exponential  factors  for  reactions  (21)-(26)  were  calculated  by 
the  procedure  described  under  COMPUTATIONAL  PROCEDURE  with  the 
parameters  of  the  reactants  in  Table  I  and  those  of  the  linear  activated  com¬ 
plexes  in  Tables  II  and  III.  The  anharmonicity  corrections  to  the  real  fre¬ 
quencies  were  neglected  since  Weston  showed  that  this  effect  is  negligible 
for  the  transition  state  when  compared  to  other  uncertainties  of  the  Sato 
method.  The  tunnelling  correction  t  was  calculated  from  the  expression 
obtained  by  Bell  (Eq.  19).  The  harmonic  force  constant  which  corresponds 
to  the  imaginary  frequency  of  the  reaction  coordinate  was  not  used  to  evaluate 
the  tunnelling  correction  factor  for  each  reaction.  Instead,  a  new  imaginary 
frequency  was  calculated  for  each  reaction  by  fitting  the  best  obtainable  para¬ 
bolic  potential  to  the  Sato  barrier  over  that  region  where  tunnelling  is  appre¬ 
ciable.  The  transmission  coefficient  of  one -half  was  used  for  isotopically 
mixed  molecules,  and  of  unity  for  all  others. 
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Figure  8  compares  the  log  obtained  experimentally  with  that 

calculated  in  this  paper.  The  theoretical  absolute  rates  for  calculated 

from  the  Sato  potential  agree  quite  well  with  that  of  the  experimental  data. 

2 

The  experimental  value  of  k, .  obtained  by  Ashmore  and  Chanmugam  at 

o  **  *  20 

523  K  compares  favorably  with  the  results  of  Rideal  and  Steiner  U  at  1000°K 

21 

and  of  Klingelhoefer  and  Rodebush  at  room  temperature.  These  three 
experimental  points  lie  on  a  straight  line;  from  the  slope  of  this  line  the 
activation  energy  of  reaction  (21)  is  5.  500  ±  0.  200  kcal/mole.  This  value  of 
the  activation  energy  compares  favorably  with  a  value  of  5.  90  kcal/mole 
calculated  by  Rideal  and  Steiner  from  their  experimental  data  at  1000°K  and 
5.  6  kcal/mole,  which  they  calculated  from  the  results  of  Klingelhoefer  and 
Rodebush. 

The  relative  rates  of  reactions  of 


21 


H  +  HCl 


H2  + 


Cl 


(21) 


and 


H  +  Cl2  “ ♦  HCl  +  Cl  (22) 

were  determined  in  the  temperature  range  273-33 5°K  by  Klein  and  Wolfsberg^. 
They  found  that  the  relative  rate  is 

^2  1 

r=-t-  =  (0.  143  ±  0.  033)  exp  -(1540  ±  130/RT)  (27) 

*22 


Klein  and  Wolfsberg  used  the  pre- exponential  factor  in  Eq.  (27)  to  determine 

the  doubly  degenerate  bending  frequency  of  the  linear  transition  state  HCl,, 

- 1  L 

and  found  it  to  have  a  value  of  about  105  ±  20  cm  .  We  evaluated  the  doubly 
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degenerate  bending  frequency  for  a  linear  transition  state  HCl^  from  the 
force  constants  (Table  III),  and  found  it  to  be  51  cm-',  The  Klein  and 
Wolfsberg  calculation  assumes  that  the  stretching  frequency  of  the  HCl^ 
transition  state  is  large  enough  to  be  neglected  in  the  temperature  range  of 
their  study  and  that  this  frequency  will  be  larger  than  the  Cl-  frequency. 

The  latter  conclusion  was  based  on  their  findings  with  respect  to  the  H^Cl 
transition  state,  where  the  most  reasonable  model  for  the  H^Cl  transition 
state  was  one  with  a  large  H-Cl  force  constant.  Our  calculation  (Table  III) 
shows  that  the  stretching  frequency  of  the  HCl^  transition  state  is  lower  than 
the  Cl.,  frequency.  In  addition,  we  calculated  the  H-Cl  and  Cl-Cl  distances 
in  the  HCl^  transition  state  to  be  r^j  =  2.  005A  and  r^  =  2.  252A. 

Klein  and  Wolfsberg  estimated  r j_j  =  1.  5A  and  r^.j  =  2.  2A.  Conse¬ 

quently  we  re-evaluated  the  bending  frequency  of  the  HCl^  transition  states, 
using  the  data  in  Tables  I  -  III  for  the  HCl^  and  H^Cl  transition  states  along 
with  the  experimental  value  of  the  pre- exponential  factor  obtained  by  Klein 
and  Wolfsberg  (Eq.  27).  The  doubly  degenerate  bending  frequency  was  found 
to  be  68  ±  20  cm  *,  which  agrees  better  with  the  value  obtained  from  the 
force  constants  calculated  from  the  Sato  potential  barrier  for  the  KCl^  linear 
transition  state.  A  comparison  of  the  relative  rates  of  reaction  {Ti)  and  (22) 
is  shown  in  Fig.  9.  The  theoretical  absolute  relative  rates  of  reaction  of 
k21^k22  ca^cu^afe<^  from  the  Sato  potential  agree  with  the  experimental  data 
within  20  percent. 

Figure  10  contains  a  plot  of  the  relative  rates  of  reaction  of  H?  and  HD 

4  c 

with  chlorine  atoms.  Bigeleisen,  et  al  ,  measured  these  rates  over  the 
temperature  range  243-350°K  and  found  that 


R  = 


21 


k25  +  k26 


(1.  24  ±  0.  03)  exp(490  ±  6/RT) 


(28) 


It  can  be  seen  from  Fig.  10  that  our  calculated  values  of  R  agree  very  well 
with  experimental  values. 
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A  plot  of  the  relative  rates  of  reaction  of  and  with  chlorine  atom 
based  on  the  parameters  given  in  Tables  I-III  is  given  in  Fig.  11.  We  note 
that  log  ^21^23  *s  *arSer  than  the  experimental  values.  Bigeleisen,  et  al, 
suggested  that  their  calculated  values  of  log  ^21^23  cou^  be  t>rought  into 
closer  agreement  with  the  experimental  values  by  choosing  the  bending 
frequency  of  the  l^Cl  complex  larger  than  200  cm'*  and  the  imaginary  fre¬ 
quency  smaller  than  720i.  Since  in  this  study  no  attempts  were  made  to  fit 
the  theoretical  structure  of  the  complex  to  the  experimental  data,  it  is  inter¬ 
esting  to  note  that  even  here  the  theoretical  values  are  within  a  factor  of  three 

3 

of  the  only  two  experimental  points  reported  by  Rollefson  .  Our  calculated 
difference  in  energies  E^-E^j  of  1.  17  1-cal/mole  compares  favorably  with 
Rollefson's  E^-E^j  of  1.  63  kcal/mole.  Rollefson  states  that  a  10-percent 
error  in  one  of  the  two  measurements  of  the  ratio  ^21^23  corresPont*s  to  a 
change  of  0.  54  kcal/mole  in  E^-E^  and  a  change  in  the  relative  rates  of  a 
factor  of  2.  7.  The  latter  would  suggest  that  additional  experiments  should  be 
carried  out  before  absolute  rate  theory  can  be  used  to  test  the  relative  rates 
of  reactions  (21)  and  (23). 

CONCLUSION 

These  theoretical  calculations  show  that  the  rates  obtained  from 
absolute  reaction  rate  theory  with  a  Sato  potential  energy  function  agree  very 
well  with  the  few  experimental  values  reported  in  the  literature.  Indeed,  they 
show  that  the  experimental  values  of  the  pre -exponential  factors  in  the 
Arrhenius -type  expression  for  rate  constants  are  not  required  to  obtain  better 
information  on  the  structure  and  vibrational  frequencies  of  the  transition 
states.  The  agreement  between  theory  and  experiment  constitutes  a  positive 
test  of  the  usefulness  of  this  transition-state  formulation  in  the  evaluation  of 
absolute  reaction  rates. 
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Figure  1.  Relative  reaction  rates  of 

isotopic  hydrogen  molecules 
with  chlorine  atoms  as  a 
function  of  temperature. 
Solid  lines  were  calculated 
by  Wheeler,  Topley,  and 
Eyring  for  the  H^Cl  tran¬ 
sition  state. 


Figure  2.  Potential  energy  surface  for  H  + 
in  vicinity  of  transition  state  usii 
Sato  method  of  construction  of  surfaces, 
K  =  0.  (Contour  lines  are  expressed 
in  kcal/mole.  ) 
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Figure  3.  Potential  energy  surface  for  H  +  HC1  in  vicinity  of  transition 
state  using  the  Sato  method  of  construction  of  surfaces, 

K  =  0.  168.  (Contour  lines  in  kcal/mole. ) 


Figure  4.  Potential  energy  surface  for  D  +  Cl,  in  vicinity  of  transition 
state  using  the  Sato  method  of  construction  of  surfaces,  K  *  0. 
(Contour  lines  in  kcal/mole. ) 


Figure  5. 


Figure  6. 


r2J ,  A 


Potential  energy  surface  for  D  +  HC1  in  vicinity  of  transition 
state  using  the  Sato  method  of  construction  of  surfaces, 

K  a  0.  168.  (Contour  lines  in  kcal/mole. ) 


Potential  energy  surface  for  H  +  DC1  in  vicinity  of  transition 
state  using  the  Sato  method  of  construction  of  surfaces, 

K  s  0.  168.  (Contour  lines  in  kcal/mole.  ) 
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Figure  7.  Potential  energy  surface  for  D,  +  Cl  in  vicinity  of  transition 
state  using  the  Sato  method  of  construction  of  surfaces, 

K  :  0.  168,  (Contour  lines  in  kcal/mole.  ) 


Figure  8.  Log^k^j  versus  10^/T 
for  reaction  H,  +  Cl 
going  to  H  +  HC1 
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Figure  9-  Relative  rates  of  reaction  of  HC1  and  Cl?  with  chlorine  atoms 


Figure  11.  Relative  rates  of  reaction  of  H,  and  D,  with  chlorine  atoms 
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